Entropically Patchy Particles 



o 

<N 

Oh 
< 

(N 



O 



o 



> 

in 

in 

o 
en 

> 

X 
S3 



Greg van Anders/ N. Khalid Ahmed/ Ross Smith/ Michael Engel/ and Sharon C. Glotzer^'^ 

^Department of Chemical Engineering, University of Michigan, Ann Arbor, MI 48109-2136, USA 
^Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI 48109-2136, USA 

Directional entropic forces that cause particle alignment have recently been proposed in systems of hard 
polyhedra that order into crystals. Here we provide a means of quantifying these forces and find them to be 
of the order of a few k^T. We show that these forces originate from "entropic patches," which are features 
of particle shape that promote local dense packing. Using the notion of entropic patches, we engineer particle 
shape to target the assembly of specific crystal structures. We show that this procedure can be generalized as 
anisotropy dimensions, similar to those exploited for enthalpically patchy particles. 



It is known since the work of Kirkwood on hard spheres (B 
and Onsager on hard rods O that entropy can order hard par- 
ticles. Recent studies have explored the role of entropy in or- 
dering hard, convex shapes between the sphere and rod limits 
|3 -20|. The role of shape in self-assembly has been reviewed 
recently in |21 1. A general observation from many of these 
studies is that hard polyhedra align their facets at sufficiently 
high packing fractions (typically, greater than 50%). Based 
on this observation, Damasceno, et al. CSlKl argued for the 
existence of "directional entropic forces" between hard parti- 
cles. Due to these forces, spontaneous first order transitions 
were observed in computer simulation where particles aligned 
upon increasing packing fraction from a disordered fluid to 
colloidal quasicrystals, liquid crystals, and crystals with unit 
cells containing as many as 52 particles [,13l[T4ll . 

Here we compute the directional entropic forces in hard 
particle systems, by introducing a potential of mean force 
and torque (PMFT). We show that the effective interaction 
strengths can be on the order of a few /cbT, putting them in 
a range that makes them an important consideration for self- 
assembly in experimental colloidal systems. We then show 
that directional entropic forces can be engineered through the 
systematic alteration of particle shape to target specific self- 
assembled structures. We conceptualize this shape engineer- 
ing as the creation of "entropic patches," in analogy with 
the engineering of enthalpic patches using chemical or other 
means introduced in |22| and reviewed in 12311241 . In con- 
trast to the charge, chemical, etc. mediation of the interaction 
between sticky patches on enthalpically patchy particles, at- 
tractive entropic patches are features in particle shape that pro- 
mote dense packing. Finally, we demonstrate the generality of 
directional entropic forces as a design concept by introducing 
anisotropy dimensions |23 1 for entropic patchiness. 

As a means of quantifying directional entropic forces, we 
examine a single pair of hard particles in a larger system of 
identical hard particles. Such a treatment was first given in 
the case of spheres by de Boer |25|; for aspherical particles, 
first steps were taken in this direction in 1 26", TT] . At constant 
volume and number of particles, the partition function of the 
system has the form 



and {Q} are the particle orientations. We have identified the 
particle pair of interest with the labels 1 and 2. All other par- 
ticles are indicated with a tilde. For axisymmetric particles 
(to which we restrict our consideration here for brevity, see 
|28 1 for the general case), the pair has only four rotational 
and translational degrees of freedom. If their centers of mass 
are separated by (fi2, and their symmetry axes are ni^2, then 
those four degrees of freedom, written to be invariant under 
rotations of the whole system, are 

{R= \q12l 01 = fii ■qi2, 02 = -^2-^12, X = -ni-n2} . 

(2) 
Fixing the values of these coordinates allows us to find F2, the 
effective (see, e.g., 1291 ) potential of mean force and torque 
(PMFT), from 

Z= /^i?d0id02^xe"^^'^^'^''^''^^ 

= / dRd(j)id(j)2dxJ{R, (1)1^2, x) 
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The Jacobian J is 
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(3) 



(4) 



Z = J dqi^2dQi,2[dq][dQ] , (1) 

where {q} are the center of mass coordinates of the particles 



H is the Heaviside step function, and d{R, 0i, 02, x) is the 
minimum separation distance of the particle pair in their rela- 
tive position and orientation, which is negative when the par- 
ticles overlap, and positive when they do not. Note that if we 
were considering particles that interacted other than sterically, 
the Heaviside function would be replaced with the Boltzmann 
factor of the pair interaction energy. F is the Helmholtz free 
energy available to other particles in the system when the rela- 
tive position and orientation of the pair is fixed. F2 is then the 
free energy of the entire system, including the particle pair of 
interest and the remaining degrees of freedom, when the rel- 
ative position and orientation of the pair are fixed. Equating 
the integrands and taking the logarithm of both sides gives the 
PMFT 

/3F2(i?,01,02,X) =m^,01,02,X) 

- \og{J{R, 01, 02, x)H{d{R, 01, 02, x)) 

(5) 





^=i ^=i 


R^O.3 
R = OA 


3.48 ±0.05 1.16 ±0.01 
3.64 ±0.06 2.39 ±0.03 



TABLE I. Angular dependence of the PMFT. Figures are height of 
the potential in k^T above the global minimum for slice of the po- 
tential with ^ = 01 = 02 = X fixed. Although the angular differ- 
ences are small (perfect alignment is V^ = 1), the effective interaction 
varies by more than 2 ksT over this angular range at small separa- 
tions, indicating that the penalty for small misalignment is signifi- 
cant. Inset particle images illustrate the relative orientations shown. 
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FIG. 1 . The PMFT exhibits strong directionality of entropic forces 
between particles, shown here through a slice of the force compo- 
nent of the PMFT with angular directions integrated out. It displays 
a force towards the tetrahedral coordination (i) and (iii) of tetrahe- 
drally facetted spheres, compared with non-tetrahedral coordination 
(ii), here shown for a fluid at a density of 50%. The dashed curve 
shows the location of the steric exclusion boundary for the center of 
mass of the second particle around the inset particle. 



The mean directional entropic forces and torques responsi- 
ble for particle alignment are encoded in the PMFT ([5]). If 
we were to label just a single pair of particles in a system of 
colloidal particles subject to thermal (Brownian) motion, we 
would infer that the force and torque on particle 1 due to par- 
ticle 2 would be given by 



dF2 dF2 

Fi = -ri2^^ - (ni - ri2^i)^- (n2 ± ^12^2) 

on o<pi 

^ . . dF2 ^. . dF2 

d(t)i dx 



dFo 



(6) 



These equations specify the forces and torques we would ex- 
pect to measure if, e.g., we grabbed a particle pair in a system 
of hard colloids using optical tweezers and attempted to pull 
them apart. 

As a simple example of the role of the PMFT in aligning 
particles, consider a particle obtained from a sphere of radius 
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FIG. 2. Self- assembled structures of facetted spheres obtained in MC 
simulation by designing particle shape. Panel (a) shows a diamond 
lattice obtained by giving a particle tetrahedrally coordinated facets 
as in (c). Panel (b) shows a simple cubic lattice, obtained by giving 
particles cubically coordinated facets as in (d). 



r by cutting away the part of the sphere that intersects the half 
space R^ for which z/r > a to create a singly facetted sphere. 
We performed MC simulations of systems of 1000 such parti- 
cles with a = 0.01 (nearly hemispherical) at fixed volume. In 
table |T| we show the role of the PMFT in generating entropic 
torques that align particle facets in this system at a density of 
50%. We find the potential difference between different parti- 
cle orientations at fixed separation distance is of the order of 
a few /cbT, giving rise to entropic torques strongly favoring 
alignment. In Fig. [T] we show the role of the PMFT in gen- 
erating directional entropic forces that cause particles to have 
positional correlation at facets. We plot the force component 
of the PMFT for tetrahedrally facetted spherical particles as a 
function of the Cartesian components of the separation vector 
of the particles in the frame of one of the particles. The plot 
shows potential minima in which the particle pairs have tetra- 
hedral coordination (i) and (iii), compared with, e.g., facet- 
edge configurations (ii); differences are on the order of a few 
/cbT. In both cases, the role of facets as entropic "patches" 
that favor facet alignment is evident. 

Entropic patches can be used to design particles for the 
purpose of self-assembling specific structures. In particular, 
that the forces are on the order of a few /cbT puts them in 
the range that is important for self-assembly. Indeed, experi- 
ments demonstrate that directional entropic forces can domi- 
nate assembly in the presence of weak enthalpic interactions 
|[T5][30l|3Tl. As an example, we consider targetting the assem- 
bly of a (tetrahedrally-coordinated) diamond lattice by tetra- 
hedrally facetting spheres; i.e. we "slice" four equal sized 
facets into the sphere at the locations of the faces of a reg- 
ular tetrahedron. If we consider a facetting amount of to 
be a perfect sphere, and 1 to be a perfect tetrahedron, then 
at a facetting amount of 0.6 the particles self-assemble a dia- 
mond lattice in MC simulation at a packing fraction of 60% 
as shown in Fig.|2^. Similarly, if we apply cubic facetting to 
a sphere by denoting a perfect sphere as and a perfect cube 
as 1, simple cubic lattices (Fig.l2b) assemble at an amount of 
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(i) patch size 



(iv) patch angle 
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(v) number of patches 



(iii) aspect ratio 
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(vi) shape composition operation J 



FIG. 3. (a) The creation and altering of entropic patches can be conceptuaUzed as anisotropy dimensions: four that are reaUzed continuously 
are (i) patch size, (ii) curvature radius, (iii) aspect ratio, and (iv) patch angle; two that are realized discretely are (v) number of patches, and 
(vi) shape composition operation, (b) Various anisotropy dimensions from panel (a) can be combined to recover familiar shapes. From top 
to bottom: combining the anisotropy dimensions of patch angle and patch number (iv and v from panel a) recovers various polyhedra, here a 
tetrahedron, a rhombicdodecahedron, and a cuboctahedron; changing the patch size (i from panel a) by combining two tetrahedra of different 
sizes through intersection yields various truncated tetrahedra; by changing the patch number (v from panel a) of a sphere that has been facetted 
with another sphere, various lock particles dimpled spheres can be realized. 



0.6 or more in MC simulations at packing fractions of 54% or 
more. Phases for other facetting amounts are shown in Fig. [5] 

The competition between the entropy available to a pair of 
particles in isolation from the rest of the system and the os- 
motic pressure of the rest of the system on the pair causing it 
to arrange more densely gives rise to the entropic ordering in 
Fig. 2. This competition is manifest in the two terms in ([5]). 
Here, — /3F = S, the entropy available to the rest of the sys- 
tem for a fixed configuration of the particle pair. Locally, S 
is an increasing function of the volume available to the rest of 
the system, meaning the remainder of the system has an en- 
tropic preference for the dense packing of the pair. This term 
is of the order of the stress tensor in the system. In the infi- 
nite pressure limit, relevant for packing, this term dominates. 
However, for assembly, which occurs at finite pressures, both 
terms contribute, and it is this competition that can lead to a 
difference between dense packing and assembled structures, 
as has been noted in several recent studies ifTTlfTsHTSlfTSl . 

Here we engineered particle shape to promote locally dense 
packing by facetting spheres, but there are many ways of al- 
tering particle shape to promote locally dense packing. The 
fact that directional entropic forces cause alignment of facets, 
thereby controlling the arrangement of particles in succes- 
sive neighbor shells, is reminiscent of the angular specificity 
of interactions conceptualized as patchy particles in ||22l |23| . 
The notion of (enthalpically) patchy particles provided a new 
means of conveying angular specificity to particle-particle 
interactions in a way that was previously only possible for 
molecules. Examples of patchy particles include Janus col- 
loids 132^361 . striped nanospheres (3711 and nanorods, I381I . 
and DNA-coated patchy particles (391, among many others 
(see the reviews ||23l[24l. However, in contrast to the chem- 
ical or other patterning that leads to enthalpic patchiness, the 
angular specificity of interactions arises here solely due to en- 







FIG. 4. (a) Many anisotropy dimensions for enthalpically patchy 
particles |23| have entropically patchy counterparts. Example tradi- 
tional sticky patchy particles left (i) and (ii) with entropically patchy 
counterparts right (i) and (ii). (b) Enthalpic patches can be combined 
with entropic patches to enhance or inhibit entropic patchiness as 
shown schematically here. This can be obtained by using the same 
direction for both types of patches (left), or by using different direc- 
tions for the entropic patches and the enthalpic patches (right). 



tropic considerations that arise from particle shape. The fea- 
tures in shape that promote dense packing are, in effect, act- 
ing as "entropic patches" that cause preferential alignment. 
We introduce a design space for these entropic patches, as 
was done for enthalpic patches (231, through anisotropy di- 
mensions (Fig. [3^). In contrast to traditional sticky patches, 
entropic patchiness is an emergent, effective concept that de- 
pends on packing fraction and on all geometric features that 
come within the characteristic length of the interaction at 
a given packing fraction. Such features may include, e.g.. 



flat facets or other low curvature regions, and interlocking or 
"mating" features. Conversely, high curvature regions could 
be used to introduce repulsive patches. 

Six examples of anisotropy dimensions are illustrated in 
Fig. [3^. Examples of continuous anisotropy dimensions in- 
clude (i) patch size, (ii) curvature radius, (iii) aspect ratio, 
and (iv) patch angle. Examples of discrete anisotropy dimen- 
sions include (v) number of patches and (vi) shape compo- 
sition operation. These anisotropy dimensions can be com- 
bined in different ways (Fig.|3|3) to recover particles of various 
known shapes, e.g. polyhedra and dimpled lock-and-key col- 
loids 1401 ED • We also note that many of these anisotropy di- 
mensions map naturally to anisotropy dimensions introduced 
for enthalpically patchy particles |23|. Examples of enthalpi- 
cally patchy particles and their entropically patchy analogues 
are shown in Fig.|4^. 

When combined with traditional enthalpic patchiness, en- 
tropic patchiness greatly enlarges the already vast design 
space for new nano- and micron- scale building blocks. The 
formalism presented here allows the assessment of the relative 
strength of entropic driving forces for assembly when other 
forces also contribute. In experimental systems, particles are 
seldom solely hard |42|. However, the fact that the scale of 
interaction is on the order a few /crT at intermediate packing 
densities means that there might be a large class of systems 
in which electrostatic or other forces can be sufficiently con- 
trolled such that entropic patches supply the dominant force 
controlling their self-assembly, as was demonstrated recently 
in ||30|l31 1. Of course, as depicted in Fig.|4]3, combining en- 
thalpic patches with entropic patches, e.g. grafting DNA or 
other ligands to facets, should enhance the tendency for face- 
to-face alignment of particles. Such interplay of forces under- 
lies bulk assemblies observed for naturally facetted nanopar- 
ticles |15|. 

Finally, we note that entropic patchiness does not require 
depletant molecules. The connection to depletion is discussed 
in a forthcoming paper |28|. 

We thank Daan Frenkel and Randall Kamien for discus- 
sions concerning the early literature. This material is based 
upon work supported by, or in part by, the U.S. Army Re- 
search Office under Grant Award No. W911NF-10-1-0518, 
the DOD/ASD(R&E) under Award No. N00244-09- 1-0062. 



Methods 

To measure the potential of mean force (PMFT) we per- 
formed Monte Carlo (MC) simulations of fluids of 1000 hard 
facetted spheres at fixed volume. Results in Table |l] are for 
nearly hemispherical particles; the radius of the sphere was 
taken to be unity, and the distance from the facet to the center 
of the sphere was taken to be 0.01 to avoid numerical errors 
that would occur in our algorithms if the center of the sphere 
was taken to lie on the the facet. The PMFT was computed 
in a fluid at equilibrium by partitioning the four scalar quan- 
tities defining independent relative positions and orientations 



into finite-size regions, and counting the relative frequency of 
observing a pair of particles in each region. Error bars are 
standard errors of the mean of independent runs of indepen- 
dently equilibrated systems. 

Similarly, for Fig. [T] MC simulations of tetrahedrally 
facetted spheres were performed. The displacement between 
particles was computed in the coordinate frame of each parti- 
cle, and possible relative positions were partitioned into finite- 
sized regions. The PMFT was computed by observing the rel- 
ative frequency of observing a pair of particles in each region. 
The results show the average of several independent runs. 

The MC method employed above, and for the results shown 
in Fig. |2] and Fig. |5] employed single particle moves for both 
translation and orientation. In all cases the simulation box was 
taken to be periodic and the volume fixed. However, for Fig. [2] 
and Fig. [5j the box was permitted to shear at fixed volume. 
In all cases, facetted spheres were approximated by polyhedra 
with 113 or more vertices, depending on the facetting amount, 
and overlaps were checked using the same implementation of 
the GJK algorithm f43l as used in fTTl. 

In Fig. [5] we show the phases that self-assembled in MC 
simulations. 
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FIG. 5. Self-assembled phases of entropically patchy hard spheres. We explore the effect of the anisotropy dimensions of facetting number, 
angle, and size from Fig.|3] Sufficiently strong entropic patches require at least some facet alignment, which in turn determines the final crystal 
structure. 



